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EVALUATION  OF  A  SIMPLE  MODEL  FOR 
PREDICTING  PHOSPHORUS  REMOVAL 
BY  SOILS  DURING  LAND  TREATMENT 
OF  WASTEWATER 


J.C.  Ryden,  J.K.  Syers  and  I.K.  Iskandar 


INTRODUCTION 

About  10  or  15  years  ago,  there  was  considerable 
controversy  over  the  relative  effect  of  external  P 
input  to  lakes  versus  internal  input  (release  ofP  from 
sediment)  on  the  excessive  fertilization  of  water 
bodies.  It  has  since  been  established,  however,  that 
reducing  the  external  P  input  to  lakes  can,  in  many 
cases  (such  as  Lake  Washington,  Seattle),  reduce 
algal  growth  in  proportion  to  the  reduction  in  the 
external  P  load.  Reducingp  concentration  in  waste- 
water  to  <0.1  ^g  P/mL  by  advanced  waste  treatment 
methods  was  necessary  to  achieve  this  reduction  in 
P  pollution. 

Land  application  of  liquid  waste,  however,  has 
recently  gained  increasing  attention  as  a  viable  al¬ 
ternative  to  stream  discharge  or  even  to  advanced 
waste  treatment  (Iskandar  et  al.  1 976).  Because  a 
major  objective  of  land  treatment  is  to  reduce  the 
concentration  of  P  in  wastewater  to  some  acceptable 
level  (perhaps  similar  to  advanced  waste  treatment), 
and  because  plant  uptake  of  p  constitutes  only  a 
small  proportion  of  the  P  applied,  a  knowledge  of 
the  soil  properties  that  facilitate  P  removal  is  re¬ 
quired.  These  factors  should  be  considered  in  the 
se'  ciion  of  a  site  for  wastewater  land  application. 
Some  soil,  such  as  that  used  in  the  Manteca,  Cali¬ 
fornia,  land  application  system,  has  a  very  limited 
capacity  to  remove  P  from  municipal  wastewater 
(Iskandar  and  Syers  1980)  and  should  have  not 
been  selected  for  such  purpose.  Soils  vary  consider¬ 
ably  in  this  ability  to  sorb  added  P,  depending  on  the 
amounts  and  nature  of  soil  components,  soil  pH, 
temperature,  water  content,  amounts  and  forms  of 
P  in  the  applied  effluent,  and  management  practices 
(Monke  et  al.  1974,  Barrow  and  Shaw  1975,  Toffie- 
mire  and  Chen  1976,  Tofflemire  et  al.  1978,  Syers 
and  Iskandar  1 981 );  sc  there  is  a  need  for  a  method 
to  characterize  the  soils  selected  for  land  application. 


Mathematical  models  can  provide  potentially  use¬ 
ful  information  on  site  selection  and  the  prediction 
of  site  longevity,  as  well  as  on  methods  for  improving 
the  management  of  land  application  systems.  Sev¬ 
eral  models  have  been  developed  with  a  view  toward 
describing  P  movement  in  soils.  These  range  in  com¬ 
plexity  and  scope  from  rather  simple,  empirical 
models  which  may  (Enfield  1974,  Ryden  and  Pratt 
1980)  or  may  not  (Taylor  and  Kunishi  1974,  Harter 
and  Foster  1976)  involve  kinetics,  to  models  based 
on  sorption  theory  and  kinetics  (Enfield  and  Bledsoe 
1975,  Enfield  and  Shew  1975,  Novak  et  al.  1975, 
Shah  etal.  1975,  Enfield  et  al.  1976,  Mansell  et  al. 
1977a,  Selim  1978).  In  some  cases,  more  complex, 
mechanistic  multiphase  models  have  been  adopted 
in  which  the  importance  of  kinetics  is  emphasized 
(Mansell  et  al.  1977b)  or  deemphasized  (Enfield  et 
al.  1977).  Some  of  these  models  (Enfield  1974, 
Enfield  et  al.  1 976)  are  not  concerned  with  P  move¬ 
ment,  per  se,  but  with  the  kinetics  of  P  sorption  by 
soils,  which  is  an  important  aspect  of  P  movement. 
Enfield  (1978)  has  reviewed  several  aspects  of  em¬ 
pirical  sorption  theory  and  mechanistic  multiphase 
models  for  P  movement  in  soils. 

One  of  the  simplest  models  proposed  to  date  is 
that  of  Harter  and  Foster  (1976).  It  uses  an  empir¬ 
ical  polynomial  sorption  equation  which  expresses 
the  relationship  between  the  amount  of  P  sorbed 
and  that  added.  It  takes  the  form: 

Y  =  A  +  BX  +  Of2  +  DX3  +  EX*  (1) 

where  Y  =  amount  of  P  sorbed 
X  =  amount  of  P  added  . 

To  obtain  this  equation,  Harter  and  Foster  (1976) 
added  increments  of  P  to  the  same  soil  sample  in  the 
laboratory.  As  suggested  by  Ryden  and  Pratt  (1980), 
however,  several  individual  but  different  additions 


of  P  may  also  be  satisfactory.  Harter  and  Foster 
found  that  polynomial  equations  gave  a  better  de¬ 
scription  of  P  sorption  isotherms  than  the  Freund- 
lich  equation.  Phosphate  breakthrough  curves  were 
then  calculated  by  using  a  simulation  model,  and  es¬ 
timates  of  site  longevity  were  obtained.  The  major 
’.dvantage  of  this  model  is  its  simplicity  in  that  it 
avoids  the  use  of  “cumbersome"  and  “inadequate” 
equations,  and  does  not  require  a  complicated  mois- 
aire  movement  computer  program  for  water  transport. 
This  approach,  however,  ignores  the  kinetics  of  P 
sorption  that  have  a  varying  effect,  depending  on  the 
soil  in  question  (Rydenetal.  1977b).  The  importance 
of  including  kinetics  in  empirical  models  has  been 
emphasized  by  Enfield  (1976),  Enfield  and  Shew 
(’975),  and  Enfield  etal.  (1976).  Furthermore,  the 
Harter  and  Foster  model  was  developed  using  only 
one  concentration  of  added  P  and,  without  the  adop¬ 
tion  of  a  generalized  theory  of  P  sorption  by  soils, 
the  model  requires  a  considerable  amount  of  labora¬ 
tory  work  to  determine  the  parameters  describing 
P  sorption.  In  addition,  there  has  been  no  evaluation 
of  the  effectiveness  of  the  model.  Nevertheless,  this 
and  similar  empirical  models  have  potential  foi  pre¬ 
dicting  P  movement  in  soil  and  require  field  evalua¬ 
tion. 

The  mechanistic  multiphase  model  developed  by 
Mansell  etal.  (1977b)  attempted  to  describe  a  much 
more  complex  system.  Chemical  kinetics  and  trans¬ 
port  theory  were  used  to  describe  the  interaction  and 
movement  of  P  in  soil.  Four  phases  of  soil  P  were 
assumed.  These  were  physically  absorbed  and  chemi¬ 
sorbed  P  as  defined  by  F  , ->  et  al.  (1977b)  and 
water-soluble  and  precipitated  phases.  The  kinetics 
of  reactions  between  any  two  of  the  four  phases  are 
considered  to  be  reversible  and  six  kinetic  reactions 
were  considered.  This  model  almost  certainly  repre¬ 
sents  the  first  multiphase,  kinetic  attempt  to  simu¬ 
late  P  movement.  The  model  was  used  by  Mansell 
et  al.  (1977b)  to  describe  the  transport  of  applied  P 
during  steady  water  flow  assuming  a  range  of  rate 
coefficients.  The  equations  used  were  as  follows: 

"Ml  =-0(k'A"+ksA)  +  p(k2B  +  k6D) 

(2) 

=k)0AN  -  {k2  +  kj)  pB  +  KApC  (3) 

^9- *3 pB-k+pC  (4) 

M  =*5  (0A)  ~k6pD  (5) 


where  6  =  volumetric  soil  water  content 
p  =  soil  bulk  density 
t  =  time 

A  =  concentration  of  P  in  solution 

8  =  amount  of  P  absorbed 

C  -  amount  of  P  immobilized  (chemisorbed) 

D  =  amount  of  P  precipitated 
N  =  constant  indicating  order  of  the  adsorp¬ 
tion  process 

k ,  =  rate  coefficient  for  adsorption 
k2  =  rate  coefficient  for  desorption 
/?3  =  rate  coefficient  for  immobilization 
=  rate  coefficient  for  mobilization 
•<f5  =  rate  coefficient  for  piecipitation 
k6  =  rate  coefficient  for  dissolution. 

If  sorption  is  considered  to  be  a  two-step  process, 
involving  a  more  physical  sorption  and  chemisorption 
(Ryden  et  al.  1976),  the  model  could  represent  a 
three-phase  system,  as  suggested  by  Enfield  (1978). 

If  chemisorbed  Prepresents  the  endpoint  of  the  in¬ 
teraction  between  P  and  soil  components,  and  pre¬ 
cipitation  is  discounted,  the  model  could  be  simplified 
to  a  two-phase  system,  but  still  involve  two  steps  for 
sorption.  As  it  stands,  the  model,  although  versatile, 
is  complex  and  requires  an  accurate  description  of 
water  transport,  and  this  is  not  always  available.  Com¬ 
bined  multiple  mechanistic  and  transport  models  may 
have  a  greater  validity  and  may  be  more  satisfactory 
in  providing  a  theoretical  description  of  P  movement 
in  soils,  but  have  yet  to  receive  adequate  field  testing 
and  evaluation. 

The  complexity  of  multiple  mechanistic  and 
transport  models,  however,  may  pose  limitations 
from  a  practical  standpoint  in  that  they  require  the 
evaluation  of  numerous  parameters  relating  to  P 
sorption  by  the  soil  at  the  treatment  site,  as  well  as 
the  measurement  of  a  range  of  field  variables.  Con¬ 
sequently,  there  is  a  need  for  a  generally  applicable 
and  simple  model  to  predict  P  movement  in  soils  that 
can  be  applied  to  existing  treatment  areas  or  areas 
under  consideration  for  use  in  land  treatment  that 
have  not  been,  or  are  unlikely  to  be,  intensively  mon¬ 
itored.  The  objective  of  the  present  study  was  to 
evaluate  a  simple  model  that  was  developed  with 
these  more  practical  limitations  in  mind. 

THE  P  MODEL 
General  concepts 

The  concept  of  the  model  to  be  evaluated  in  the 
present  study  has  been  discussed  in  detoil  by  Ryden 
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Figure  I.  Generalized  equilibrium  isotherm  for  P  sorption  by  contrasting 
soils,  where  0  (the  fractional  saturation  of  the  equilibrium  sorption  maximum) 
is  plotted  against  the  equilibrium  solution  P  concentration.  The  data  points 
indicated  are  taken  from  Ryden  etal.  (1976)  and  Ryden  and  Pra  tt  (1 980). 

The  solid  line  delineating  the  isotherm  is  that  described  by  eq  8.  (•)  Ramona 
soil;  (m)  Waikakahi  soil;  (a)  Egmont  soil;  (*)  Porirua  soil;  and  (A)  Okaihau  soil. 


and  Pratt  (1980).  The  model  relies  on  the  conform¬ 
ity  of  P  sorption  by  soil  at  the  treatment  area  to  a 
generalized  isotherm  (Fig.  1)  describing  Psorption 
by  soils  (Ryden  and  Syers  1977a).  Conformity  of 
P  sorption  to  this  generalized  isotherm  has  been 
observed  for  a  wide  range  of  contrasting  subsoils 
and  topsoils  (Ryden  et  al.  1976,  Ryden  and  Syers 
1977a,  Ryden  and  Pratt  1980).  The  generalized 
isotherm  describes  the  relationship  between  the 
fractional  saturation  ( 0 )  of  the  maximum  Psorption 
ability  of  the  soil  and  the  solution  P  concentration 
at  equilibrium.  This  isotherm  provides  a  simple 
method  of  assessing  equilibrium  sorption  capacity 
at  a  specified  solution  P  concentration.  The  specified 
solution  P  concentration  may  be  as  high  as  that  in 
the  wastewater  applied  or  as  low  as  that  regarded  as 
being  a  maximum  by  a  regulatory  agency.  Thus  the 
intention  of  the  model  is  to  calculate  the  capacity 
of  a  soil  to  sorb  P  at  any  given  equilibrium  solution 
P  concentration. 

From  this  estimate  of  P  sorption  capacity  at  the 
specified  solution  p  concentration  (in  most  cases, 
the  mean  dissolved  inorganic  P  concentration  of  the 
applied  effluent),  the  longevity  of  the  treatment  area 
may  be  calculated  from  eq  6 


when  T  is  the  time  in  years  for  the  P  front  to  reach 
a  giver  depth  in  the  profile,  5p  is  the  P  storage  capac¬ 
ity  in  kg  ha-1  of  the  volume  of  soil  above  that  depth 
for  the  specified  solution  P  concentration,  /p  is  the  P 
input  in  kg  ha'1  yr*1  and  is  the  amount  of  P  re¬ 
moved  in  any  harvested  crop. 

Equation  6  assumes  that  the  net  P  input  (/p  -  Hp ) 
reacts  progressively  with  successive  depth  elements 
of  the  soil  profile.  The  sorption  capacity,  as  defined 
above,  of  the  first  depth  element  becomes  saturated 
before  Appreciable  movement  of  Pinto  the  next 
depth  element.  By  implication,  therefore,  the  boun¬ 
dary  between  P-enriched  and  nonenriched  soil  is 
assumed  to  be  rather  abrupt.  The  model  also  con¬ 
siders  that  water  movement  is  unimportant  relative 
to  the  kinetics  of  the  Preactions  and  that  there  is 
sufficient  time  for  slow  reactions  to  have  an  apprec¬ 
iable  impact  on  the  extent  of  Pretention  within  the 
soil  profile. 

Procedures  for  obtaining  model  parameters 

Values  for  the  parameters  in  eq  6  are  relatively 
easy  to  obtain.  The  values  of  lp  and  Hp  can  be 
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determined  with  reasonable  accuracy  by  well  estab¬ 
lished  techniques.  The  value  of  Sp  is  estimated  from 
the  generalized  P  sorption  isotherm  (Fig.  I),  as  de¬ 
scribed  below. 

The  procedures  developed  for  use  in  the  present 
study  involve  the  determination  of  the  amount  of 
P  sorbed  at  equilibrium  for  additions  of  P,  resulting 
in  equilibrium  P  concentrations  of  between  1  and 
10  mg  PL'1.  The  amounts  of  P  to  be  added  should 
be  no  more  than  600,  1 500,  and  3000  pg  P  g'1  soil 
for  soils  with  expected  “low,”  "medium,"  and  "high” 

P  sorption  capacity,  respectively.  A  qualitative  esti¬ 
mate  of  sorption  capacity  can  be  made  from  available 
information  on  the  mineralogical  composition  of  the 
soil,  and  to  some  extent  from  soil  color,  which  for  B 
horizon  material  should  provide  an  indication  of  the 
amount  of  hydrous  iron  oxides.  The  laboratory  studies 
to  determine  the  value  of  Sp  must  also  use  soils  that 
have  received  no  wastewater.  These  can  be  collected 
from  an  area  prior  to  wastewater  treatment  or  from 
zones  contiguous  with  an  existing  treatment  area. 

From  the  determinations  of  equilibrium  P  sorption, 
conformity  to  the  generalized  isotherm  can  be  assessed, 
and  the  P  sorption  capacity  at  a  specified  solution  P 
concentration  can  be  calculated.  The  precise  proce¬ 
dure  is  described  below. 

Samples  (4  g)  of  air-dried,  <2-mm  soil  are  shaken 
with  400  mL  of  10'-  M  CaCI2  containing  ^0  mg 
L'1  Hg Cl2  and  three  to  five  different  levels  of  added 
P  as  KH2P04.  Shake  the  bottles  on  an  end-ovc'-end 
shaker  for  approximately  72  hr.  Remove  the  bottles 
and  allow  to  settle  for  about  2  hr  or  until  the  super¬ 
natant  solution  is  clear.  Remove  a  4-mL  aliquot  of 
the  dear  supernatant  from  the  upper  2  cm  for  the 
determination  of  inorganic  P  (Murphy  and  Riley 
1962).  Recap  the  bottles,  continue  shaking,  and 
repeat  the  sampling  procedure  at  approximately  98, 

1 20,  144,  168  and  240  hr  using  the  procedure  de¬ 
scribed  above.  Solution  P  concentration  at  each 
time  interval  is  plotted  against  the  reciprocal  of  time, 
and  what  should  be  a  linear  relationship  is  extrapo¬ 
lated  to  infinite  time,  t  =  °°,1  It  =  0  (Ryden  and 
Syers  1975).  The  solution  P  concentration  at  /  =  00 
is  taken  as  the  equilibrium  value.  From  this  value 
and  the  values  for  the  amount  of  P  added  and  the 
soil-to-solution  ratio,  the  corresponding  value  for 
“equilibrium  sorbed  P”  is  calculated. 

The  amount  of  native  sorbed  P  present  in  the  soil 
before  P  addition  must  be  added  to  the  value  of  equil¬ 
ibrium  sorbed  P  to  give  an  estimate  of  total  equilibrium 
sorbed  P  (xe)  (Ryden  etal.  1977b).  The  amount  of 
native  sorbed  P  is  determined  by  shaking  the  air-dried, 
<2-mm  soil  (2  g)  with  30  mL  of  0.1  M  NaOH  for  20 
hr.  The  suspension  is  centrifuged  and  the  supernatant 
decanted  off.  Any  soluble  organic  material  in  solution 
is  removed  using  P-free  activated  charcoal.  A  suitable 
aliquot  is  neutralized,  and  dissolved  inorganic  P  is  de¬ 


termined.  The  amount  of  P  extracted  is  expressed 
per  unit  weight  of  soil. 

The  equilibrium  sorption  maximum  (be)  may  then 
be  calculated  (eq  7)  from  xt  and  by  reference  to  the 
generalized  isotherm  to  obtain  the  value  of  d  for  the 
observed  equilibrium  solution  P  concentration 

be=X'!0  (7) 

where  0  is  the  fractional  saturation  of  the  sorption 
maximum  for  the  experimentally  determined  equil¬ 
ibrium  solution  P  concentration.  The  value  of  6  can 
be  most  accurately  determined  from  eq  8  which  de¬ 
scribes  the  generalized  isotherm  (Fig.  1 )  over  the  solu¬ 
tion  P  concentration  ( c )  range  of  0  to  1 2  mg  L'1 . 

0  =  0.778  -  (0.455e'°-148‘) 

-  (0.217e-7-56c)  .  (8) 

Conformity  of  P  sorption  by  the  soil  under  evalu¬ 
ation  to  the  generalized  sorption  isotherm  can  be 
verified  by  comparison  of  the  £>e  values  calculated 
from  eq  7  for  the  different  values  of  xe  obtained 
during  the  sorption  study.  Conformity  to  the  gen¬ 
eralized  isotherm  is  verified  by  essentially  constant 
be  values. 

When  conformity  to  the  generalized  isotherm  has 
been  demonstrated  for  one  depth  element  within  the 
soil  profile,  i!  is  reasonable  to  assume  that  the  re¬ 
maining  depth  elements  will  also  conform  to  this 
sorption  model.  Subsequent  determinations  of  be 
for  other  depth  elements  will  therefore  require  only 
one  level  of  P  addition.  This  assumption,  however, 
must  not  be  applied  if  the  soil  profile  under  consider¬ 
ation  shows  a  discontinuity  in  its  mineralogical  compo¬ 
sition. 

The  value  of  xe  corresponding  to  the  specified  so¬ 
lution  P  concentration  for  each  depth  element  may 
then  be  determined  from  eq  7  using  the  estimate  of 
bt,  obtained  experimentally  as  described  above,  and 
the  value  of  6  calculated  from  eq  8  for  the  specified 
solution  P  concentration.  The  value  xt  minus  the 
amount  of  native  sorbed  P  may  then  be  set  asS  when 
the  extent  of  P  sorption  is  expressed  on  a  weight/vol¬ 
ume  basis,  which  requires  knowledge  of  the  bulk  den¬ 
sity  of  the  soil  under  investigation. 

EVALUATION  OF  THE  MODEL 

The  submodel  describing  P  sorption  has  been 
tested  for  a  wide  variety  of  soils  with  various  objec¬ 
tives  in  view  (Ryden  et  al.  1976,  Ryden  and  Syers 
1977a,  Ryden  and  Pratt  1980).  The  use  of  this  sub¬ 
model  within  the  P-balance  model  described  by  eq 
6  has  not  been  previously  evaluated.  An  application 
that  attempts  to  predict  storage  capacity  and  site 
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longevity  at  two  different  sites  in  New  England  with 
three  different  soils  is  presented  below. 

Description  of  study  sites 

The  land  treatment  sites  used  for  model  evalua¬ 
tion  in  the  present  study  are  located  at  the  Cold 
Regions  Research  and  Engineering  Laboratory 
(CRREL)  in  Hanover,  New  Hampshire,  and  at  Fort 
Devens,  Massachusetts.  The  CRREL  facility  consists 
of  six  outdoor  test  cells  (8.4x  8.4  m  and  1 .5  m  deep) 
containing  either  Windsor  sandy  loam  or  Charlton 
silt  loam  soils.  Soil  horizons  are  0  to  1  5  cm  (A), 

1 5  to  45  cm  (B),  and  45  to  1 50  cm  (C)  for  both 
soils.  Particle  size  analysis  of  these  soils  has  been 
presented  elsewhere  (Iskandar  et  al.  1979).  Waste- 
water  applied  to  the  test  cells  containing  Windsor 
and  Charlton  soils  had  received  preapplication  treat¬ 
ment  to  produce  either  primary  or  secondary  efflu¬ 
ents  followed  by  on-line  disinfection  with  ozone  or 
chlorine.  Wastewater  application  commenced  in  the 
spring  of  1973  and  was  applied  by  spray  irrigation 
in  varying  quantities  during  the  growing  season. 

During  the  winter  of  1974/1975,  wastewater  was 
also  applied  by  inundation.  A  mixture  of  forage 
grass  was  planted  in  1973;  however,  quackgrass 
(Agropyron  repens  L.)  was  the  predominant  species 
during  1974  and  1975.  Forages  were  harvested 
three  to  four  times  each  year.  Amounts  of  P  applied, 
effluent  P  concentration  and  plant  uptake  of  P  be¬ 
tween  1973  and  1978  are  discussed  later  (see  Table  3). 

For  model  evaluation,  samples  of  the  soil  solution 
were  removed  using  suction  lysimeters  installed  in 
one  test  cell  for  each  soil  at  mean  depths  of  9.8,  40 
and  75  cm.  Soil  samples  were  also  collected  from 
each  test  cell  during  1979  and  the  amount  of  sorbed 
P  was  determined  by  extraction  with  0.1  M  NaOH 
as  described  above. 

The  Fort  Devens  site  is  a  rapid  infiltration  treat¬ 
ment  facility  which  has  been  in  operation  since  1946 
serving  a  U.S.  military  installation.  The  system  con¬ 
sists  of  22  treatment  beds.  In  1973,  the  flow  rate 
ranged  from  2676  to  9541  m3  day"1  (1.3  mgd). 

Every  ?4  days,  15  to  20  cm  of  primary-treated,  un¬ 
disinfected  wastewater  is  applied  to  each  infiltration 
basin.  The  wastewater  infiltrates  within  2  to  3  days. 
During  the  summer,  vegetation  consists  of  weedy 
grasses  which  include  (Satterwhite  et  al.  1976)  fall 
panicum  (panicum  dicotomilflorium  Michx.)  and 
barnyard  grass  (Echinochloa  crusgaUi  (L.)  Beauv). 

The  soil  at  this  treatment  site  is  a  gravelly  sand  with 
interspersed  lenses  of  silty  sand  and  sandy  gravels. 
Particle  size  analysis  indicated  that  only  1  to  3%  of 
the  particles  are  <0.054  mm,  40%  are  gravel  and  the 
rest  sand  (Satterwhite  et  al.  1976).  The  amounts  of 
P  applied  and  effluent  P  concentration  are  discussed 


later  (Table  3).  Soil  samples  were  collected  from  the 
treatment  site  during  1979.  The  amounts  of  sorbed 
P  were  determined  by  extraction  of  P  with  0.1  M 
NaOH  for  several  depth  elements  within  the  profile, 
down  to  approximately  1 .6  m. 

Conformity  of  P  sorption  to  generalized  isotherm 

The  time  dependence  of  P  sorption  by  surface 
horizons  of  the  Windsor,  Charlton,  and  Fort  Devens 
control  soils  during  periods  in  excess  of  72  hr  is 
illustrated  in  Figure  2.  For  each  soil  and  P  addition, 
an  essentially  linear  relationship  was  obtained  which 
was  extrapolated  to  infinite  time  (1  //  =  0)  to  provide 
an  estimate  of  the  equilibrium  solution  P  concentra¬ 
tion.  The  equilibrium  solution  P  concentration 
varied  widely  between  soils  and  with  the  level  of  P 
addition,  but  were  all  within  the  desired  range  of 
approximately  1  to  10  mg  L"^ . 

The  equilibrium  solution  P  concentration  obtained 
for  each  of  four  P  additions  to  the  surface  horizons 
of  each  soil  are  given  in  Table  1 .  The  values  for 
equilibrium  P  sorption  were  calculated  from  the 
equilibrium  solution  P  concentration,  the  level  of 
P  addition,  and  the  amount  of  native  sorbed  P  present 
in  each  control  soil.  The  values  obtained  are  those 
for  xe  in  eq  7.  The  values  of  0  for  each  equilibrium 
solution  P  concentration  were  calculated  from  eq  8 
describing  the  generalized  P  sorption  isotherm  (Fig.  1 ) 
and  the  corresponding  bt  values  obtained  from  eq  7 
are  given  in  Table  1. 

Despite  the  wide  range  in  P  addition  to  each  soil, 
the  values  for  be  were  essentially  constant  for  a  par¬ 
ticular  soil.  The  coefficients  of  variation  associated 
with  the  mean  value  of  be  for  each  soil  ranged  from 
only  +2.9  to  8.5%.  These  findings  are  consistent  with 
those  for  a  Ramona  sandy  loam  from  California 
(Ryden  and  Pratt  1980).  The  consistency  of  bt 
values  for  a  particular  soil  (Table  1 )  clearly  indicates 
that  P  sorption  by  each  soil  conformed  to  the  gener¬ 
alized  isotherm  for  P  sorption  by  soils  (Fig.  1 ).  It  is 
reasonable  to  assume,  therefore,  that  the  isotherm 
can  be  used  in  the  prediction  of  5p  for  these  and 
other  depth  elements  of  the  soil  profiles  used  in  the 
present  study. 

Estimates  of  be  for  additional  depth  elements  of 
the  Windsor,  Charlton,  and  Fort  Devens  profiles  are 
given  in  Table  2.  These  estimates  were  based  on  the 
determination  of  an  xe  value  from  a  single  addition 
of  P  to  each  soil.  Each  P  addition  resulted  in  an 
essentially  linear  relationship  between  solution  P 
concentration  and  the  reciprocal  of  time.  The  values 
of  bt  for  each  depth  element  suggest  appreciable 
variation  in  equilibrium  P  sorption  within  each  pro¬ 
file,  with  highest  values  generally  being  observed  in 
the  upper  part  of  the  profile.  In  the  case  of  the 
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Table  1 .  Data  for  equilibrium  sorption  of  P  by  the  surface  horizons 
of  Windsor,  Charlton,  and  Fort  Devens  soils  derived  from  the  kinetics 
of  P  sorption  and  the  generalized  P  sorption  isotherm  for  soils. 


A  dded  P 
(ugg'1) 

Equilibrium  P 
concentration 
(mg  L~’) 

Equilibrium  P 
'  sorpiionf 

(pgg~!) 

e ** 

Equilibrium  sorption 
maximum  (bf) 

(ugg~’> 

Windsor  sandy  loam 

625 

2.10 

708 

0.559 

1270 

875 

3.45 

823 

0.641 

1280 

1250 

6.12 

931 

0.723 

1290 

1750 

9.80 

1060 

0.763 

1390 

1310i  60 

Charlton  silt  loam 

625 

1.80 

720 

0.535 

1350 

875 

2.90 

860 

0.612 

1410 

1250 

5.55 

975 

0.712 

1370 

1750 

9.35 

1090 

0.760 

1430 

1390  i  40 

Fort  Devens  gravelly  sand 

250 

0.75 

268 

0.427 

628 

750 

4.20 

423 

0.670 

631 

1125 

6.85 

S33 

0.736 

724 

1500 

10.3 

563 

0.765 

735 

680  1  58 

*  Estimated  from  Figure  2. 

*  Includes  native  sorbed  P  (293,  275  and  93  tig  P  g"1  for  Windsor,  Charlton 
and  Fort  Oevens,  respectively). 

•’Calculated  from  eq  8  describing  the  generalized  P  sorption  isotherm 
(Fig.  D. 


Table  2.  Estimation  of  bt  values  for  the  additional  horizons  of  the  Windsor, 
Charlton,  and  Fort  Devens  soils. 


Soil 

P  added 
(ug  g-] ) 

Equilibrium 
solution  P 

concentration 
(mg  L-l) 

P  sorbed 
(ggg'1 ) 

Votive 
sorbed  P 
(ggg'1 ) 

*e 

(Mg-1/ 

e 

bc 

(pgg-’l 

Windsor  B 

1000 

4.15 

585 

223 

808 

0.670 

1200 

Windsor  C 

500 

2.15 

285 

124 

409 

0.563 

726 

Charlton  B 

2000 

2.00 

1800 

142 

1940 

0.551 

3S20 

Charlton  C 

Fort  Devens 

1000 

3.45 

655 

206 

861 

0.641 

1340 

(6-25  cm) 

Fort  Devens 

750 

2.05 

495 

123 

618 

0.555 

1110 

(50-75  cm) 
Fort  Devens 

750 

4.35 

315 

52 

367 

0.678 

541 

(102-163  cm) 

750 

4.85 

265 

54 

319 

0.693 

460 

F igure  2.  Relationships  between  solution  P  concentration  and  recip¬ 
rocal  of  time  for  various  additions  of  P  (indicated  against  each  rela¬ 
tionship)  to  soil  from  the  surface  horizons  of  Windsor  sandy  loam 
(*),  Charlton  silt  loani(m),  and  Fort  Devens  gravelly  sand  (w).  Extrap¬ 
olation  (dashed  line)  to  infinite  time  (l/l  =  0)  gives  an  estimate  of 
equilibrium  solution  P  concentration. 


Charlton  and  Fort  Devens  soils,  the  second  depth 
element  has  an  appreciably  higher  P  sorption  capac¬ 
ity  than  any  of  the  other  depth  elements. 

Determination  of  5p  and  site  longevity 

As  indicated  previously,  the  value  of  5p  (eq  6) 
may  be  determined  from  bt  and  the  value  of  6  for  a 
specified  solution  P  concentration.  In  the  present 
evaluation,  the  specified  solution  P  concentration 
was  set  as  that  in  the  effluent  being  applied  to  the 
soil  profile.  These  values  are  given  in  Table  3.  Hence, 
the  storage  capacity  of  the  profile  is  assumed  to  be 
limited  by  the  xt  value  corresponding  to  the  average 
solution  P  concentration  in  the  effluent. 

The  values  of  the  parameters  involved  in  the  de¬ 
termination  of  Sp  for  each  depth  element  in  each 


profile  are  given  in  Table  3  in  the  sequence  in  which 
they  are  required  in  the  calculation.  In  the  case  of 
Windsor  soil,  the  mean  P  concentration  in  the  applied 
effluent  is  5.2  mg  L"1 .  This  concentration  gives  a  6 
value  of  0.704  from  eq  8  and  an  xt  value  of  921  pg 
P  g"^  for  the  A  horizon.  The  net  sorption  capacity 
(net  xe)  at  5.2  mg  L"1  is  628  pg  P  g-1  when  the 
amount  of  native  sorbed  P  is  deducted  from  the  value 
of  xt.  Based  on  a  depth  element  of  1 5  cm  and  a 
bulk  density  within  this  element  of  1 .38  g  cm"3,  the 
storage  capacity  Sp  is  found  to  be  1 300  kg  P  ha"1 . 
Similar  calculations  provided  estimates  of  5p  for  the 
other  depth  elements  of  each  soil  (Table  3). 

The  tout  storage  capacity  of  each  profile  can  be 
determined  by  summation  of  the  individual  Sp  values 
for  each  depth.  The  total  storage  capacity  for  the 
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Table  3.  Annual  input  and  crop  removal  of  P  and  data  required  in  the  calculation  of  stor¬ 
age  capacity  (Sp)  for  each  soil. 


Horizon 

Depth 

(cm) 

0  for  P 

concentration 

of  effluent 
(m  g-’) 

b, 

(m  r1) 

V 

(gg  g'1) 

Native 
sorbed  P 
fw  9~’) 

Net  xe 

(m  r1 ) 

Bulk 

density 

(g/cm3) 

Sp 

(kg  PI ha) 

Windsor  sandy  loam 

* 

A 

0-15 

0.704 

1310 

921 

293 

628 

1.38 

1300 

B 

15-46 

0.704 

1210 

848 

233 

626 

1.58 

2970 

C 

46-137 

0.704 

726 

511 

124 

387 

1.5  3 

6220 

Charlton  silt  loamf 

A 

0-15 

0.704 

1390 

979 

275 

704 

1.08 

1 140 

8 

1  ."—46 

0.704 

3520 

2480 

142 

2340 

1.46 

10, .,00 

C 

46-137 

0.704 

1330 

946 

206 

740 

1.71 

13,300 

Fort  Devens  gravelly  sand** 

- 

0-6 

0.764 

680 

520 

93 

427 

1.4ft 

169 

- 

6-25 

0.764 

1110 

851 

123 

728 

1.4 

910 

- 

50-76 

0.764 

541 

413 

52 

361 

1.4 

617 

* 

102-163 

0.764 

460 

351 

54 

297 

1.4 

1190 

•Average  P  cone,  of  effluem  =  5.2  mg  L"' ;  total  P  input  =  243  kg  ha"'  yr*' ;  P  removed  in  harvested  crop 
=  33  kg  ha"'  yr"' . 

t  Average  P  cone,  of  effluent  =  5.2  mg  L"' ;  total  P  input  =  145  kg  ha"'  yr*' ;  P  removed  in  harvested 
crop  =  29  kg  ha"'  yr"'. 

••Average  P  cone,  of  effluent  =  10  mg  L"' ;  P  input  =  2470  kg  ha"'  yr"' ;  P  removed  in  harvested  crop 
=  0  kg  ha"'  yr"'. 

tfEstimated  value  for  <2  mm  material  which  is  47%  of  a  given  soil  volume. 


Charlton  soil  to  a  depth  of  1 37  cm  (24,700  kg  ha"' ) 
greatly  exceeded  that  (1 0,500  kg  ha"' )  for  the  Wind¬ 
sor  soil.  This  reflects  the  appreciably  higher  sorption 
maxima  of  the  B  and  C  horizons  in  the  Charlton  soil. 

The  longevity  of  each  profile,  before  appreciable 
P  breakthrough  at  the  maximum  depth  sampled,  can 
be  calculated  from  eq  6  and  the  information  in  Table 
3.  In  the  case  of  the  Windsor  profile,  5p  =  1 0,500 
kg  ha"'  and  (/p  -  Hp)  =  210  kg  ha"'  yr"' .  Assuming 
the  continuance  of  present  management,  the  front 
of  P  enrichment  within  the  profile  should  not  reach 
137  cm  for  approximately  50  years  from  the  incep¬ 
tion  (1974)  of  land  treatment  of  wastewater  on  this 
soil.  The  longevity  for  the  Charlton  soil  to  the  same 
depth  is  213  years,  reflecting  the  appreciably  greater 
storage  capacity  of  this  profile. 

Similar  calculations  for  the  Fort  Devens  profile 
indicate  a  net  storage  capacity  to  163  cm  of  5280 
kg  ha"',  when  it  is  assumed  thatSp  values  for  depth 
elements  that  were  not  sampled  are  the  average  of 
the  values  for  the  elements  immediately  above  and 
below.  Hence,  the  longevity  of  the  Fort  Devens  pro¬ 
file  for  a  P  input  of  2470  kg  ha"'  yr"1  is  approxi¬ 
mately  2  years.  As  wastewater  has  been  added  to 
this  profile  since  1946,  these  calculations  suggest 


that  the  storage  capacity  or  the  profile  to  be  depth 
of  163  cm  has  long  been  saturated.  Field  measure¬ 
ments  of  the  P  enrichment  within  the  profile  at  the 
Fort  Devens  treatment  area  are  consistent  with  this 
prediction,  as  discussed  below. 

Prediction  of  the  depth  of  P  enrichment 

The  accumulation  of  NaOH-extractable  P  within 
the  soil  profiles  receiving  wastewater  may  be  used 
as  an  index  of  the  depth  to  which  the  front  of  P  en¬ 
richment  has  moved.  Tr.is  assumption  is  reasonable 
since  the  major  portion  of  the  total  P  applied  is  in 
the  form  of  dissolved  inorganic  P  which  will  be  largely 
retained  in  the  pool  of  sorbed  P.  Mineralization  of 
applied  organic  P  will  also  ultimately  contribute  to 
the  pool  of  sorbed  P.  It  is  this  pool  which  is  removed 
during  NaOH  extraction  (Williams  et  al.  1967,  Ryden 
et  al.  1977a). 

In  the  case  of  the  Windsor  soil,  appreciable  enrich¬ 
ment  of  NaOH-P  had  occurred  to  a  depth  of  1 5  cm 
with  additional  limited  enrichment  in  the  15  to  30 
cm  depth  (Table  4).  The  total  NaOH-P  enrichment 
of  the  profile  was  1260  kg  ha"'  which  compared  fa¬ 
vorably  with  the  net  P  input  of  1 050  kg  ha"'  be¬ 
tween  1974  and  1979.  The  data  in  Table  4  suggest 
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Table  4.  Enrichment  in  NaOH-extractable  P  in 
different  depth  elements  of  the  Windsor,  Charl¬ 
ton  and  Fort  Devens  profiles  after  treatment 
with  wastewater,  and  the  corresponding  esti¬ 
mates  of  storage  capacity. 


NaOH-extractable  P  (ky  ha~* ) 

Depth  Storage 


(cm)  Control”  1979 _ •  -  NaOH-P  capacity  t 


Windsor  sandy  loam 

0-7.5 

303 

830 

527 

650 

7.5-15 

303 

845 

542 

650 

15-30 

528 

710 

182 

1480 

30-45 

528 

517 

- 

1480 

45-60 

287 

298 

1  1 

1010 

V  _ 

1260 

Net  P  input  = 

1050 

Charlton  silt  loam 

0-7.5 

222 

956 

734 

570 

7.5-15 

222 

228 

- 

570 

15-30 

330 

428 

98 

5120 

30-45 

330 

367 

37 

5120 

1  = 

859 

Net  P  input  = 

580 

Fort  Devens  gravelly  sand 

0-6 

35 

142 

107 

169 

6-25 

150 

745 

605 

910 

50-76 

78 

680 

602 

617 

102-163 

190 

1280 

1090 

1190 

*  Soil  materials  collected  prior  to  wastewater  addition . 
t  5p  value  for  the  solution  P  conce  itrations  specified  in 
Table  3. 


that  the  front  of  P  enrichment  is  at  about  1  5  cm, 
with  some  spread  into  the  1  5-  to  30-cm  depth. 
Significantly,  however,  85%  of  the  total  P  enrich¬ 
ment  occurred  in  the  0-  to  1 5-cm  depth.  If  the  zone 
of  P  enrichment  was  defined  by  a  sharp  boundary, 
as  is  assumed  in  the  P  budget  model  described  by 
eq  6,  the  front  of  P  enrichment  would  have  been 
predicted  to  have  reached  1 2  cm,  based  on  the  net 
P  input  and  the  storage  capacity  of  the  profile.  Con¬ 
sequently,  the  model  suggests  close  agreement  be¬ 
tween  predicted  and  measured  depths  of  P  enrich¬ 
ment  in  the  Windsor  soil. 

Similarly,  the  predicted  and  measured  depths  of 
P  enrichment  were  in  close  agreement  for  the  Charl¬ 
ton  soil.  Appreciable  accumulation  of  NaOH-extract- 
abl"  P  was  only  observed  for  the  0  to  7.5  cm  pro¬ 
file  depth  (Table  4)  in  which  86%  of  the  total  P  en¬ 
richment  occurred.  Based  on  the  storage  capacity 
of  the  profile  and  the  net  P  input,  the  predicted 
depth  of  the  front  of  P  enrichment  is  7.6  cm,  in 
very  -Jose  agreement  with  the  field  observations. 


Table  5.  Measured  and  predicted  soil  solution  P 
concentrations  in  different  depth  elements  of 
the  Windsor,  Charlton  and  Fort  Devens  soils. 


Depth 

(cm) 

Dissolved  inorganic  P  in  soil  solution  (mu  L'!  1 

Measured  values 
( mean  and  range)* 

Predicted  values  f 

Windsor  sandy  loam 

7.5-12 

0.245  (0.040-0.441) 

3.0 

35-45 

0.031  (0.001-0.069) 

0.055 

70-80 

0.015  (0.001-0.060) 

n.d.** 

Charlton  silt  loam 

7.5-12 

0.041  (0  00 1  -0.067) 

0.060 

35-45 

0.019  (0.001-0.052) 

0.010 

70-80 

0.015  (0.001-0.029) 

n.d. 

Fort  Devens  gravelly  sand 

6-25 

8  (5-13) 

1.8 

50-76 

8  (5-11) 

6.8 

102-163 

8  (5-12) 

5.1 

*  On  nine  days  during  October  and  November  1979  (total 
of  27  samples),  except  Fort  Devens  for  which  values  are 
drawn  from  Satterwhite  et  al.  (1976). 
t  Predicted  from  the  existing  saturation  of  the  sorption 
maximum  after  application  of  wastewater;  details  in  text. 
**Not  determined;  appropriate  information  not  available. 


The  depth  of  P  enrichment  is  also  reflected  in 
measured  dissolved  inorganic  P  concentrations  in 
the  soil  solution  at  various  depths  in  the  undisturbed 
profile  (Table  5).  In  the  case  of  Windsor  soil,  the 
soil  solution  P  concentration  at  a  mean  depth  of  9.8 
cm  (within  the  depth  of  P  enrichment)  was  approxi¬ 
mately  10  times  greater  than  that  at  40  cm,  a  depth 
outside  the  zone  of  P  enrichment.  For  Charlton 
soil,  the  soil  solution  P  concentration  at  a  mean  depth 
of  9.8  cm  was  only  about  twice  that  observed  at  lower 
depths  where  no  P  enrichment  had  occurred.  This 
is  consistent  with  the  fact  that  the  predicted  and 
measured  depths  of  P  enrichment  were  at  only  7.5  cm 
in  the  Charlton  soil,  and  suggests  a  sharp  boundary 
between  the  P-enriched  and  nonenriched  zones  within 
the  profile. 

The  P  sorption  submodel  and  the  values  for  NaOH- 
cxtractable  P  can  also  be  used  to  predict  the  dissolved 
inorganic  P  concentrations  in  the  soil  solution  at  a 
particular  depth.  This  is  achieved  by  setting  8  = 
NaOH-P/6e  and  calculating  the  corresponding  solution 
P  concentration  from  eq  8.  Where  comparisons  can 
be  made,  predicted  solution  P  concentrations  agreed, 
with  one  exception,  moderately  well  with  the  meas¬ 
ured  values  (Table  5). 

In  the  case  of  Fort  Devens  soil,  NaOH-extractable 
P  showed  considerable  enrichment  for  all  depth 
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elements  as  compared  to  the  control  soils  (Table  4). 

The  net  enrichment  of  NaOH-P  agreed  fairly  well 
with  the  storage  capacity  of  each  depth  element  of 
this  profile,  particularly  at  the  greater  depths.  These 
data  (Table  4)  suggest  that  the  Fort  Devens  profile 
is  saturated,  with  respect  to  P  sorption  from  effluent 
having  a  mean  P  concentration  of  10  mg  L"1 .  The 
saturation  of  the  Fort  Devens  profile  was  predicted 
in  the  calculation  of  site  longevity,  and  is  reflected 
by  the  mean  and  range  of  soil  solution  P  concentra¬ 
tions  measured  at  various  depths  in  the  profile  (Table  5). 
These  concentrations  are  essentially  the  same  as  that 
for  the  untreated  effluent  and  suggest  the  lack  of 
further  P  removal  in  the  upper  163  cm  of  the  profile. 

As  for  Windsor  and  Charlton  soils,  the  soil  solution 
P  concentration  predicted  from  the  P  sorption  sub¬ 
model  were  of  the  same  order  as  the  measured  con¬ 
centrations. 


ADVANTAGES  AND  LIMITATIONS 
OF  THE  MODEL 

The  objective  of  the  model  described  and  evalu¬ 
ated  in  the  preceding  disucssion  was  to  predict  the 
storage  capacity  of  a  soil  profile  with  respect  to  P 
sorption,  and  thereby  to  predict  site  longevity  during 
land  treatment  of  wastewaters.  In  this  respect,  the 
model  has  been  successful  in  predicting  the  lower 
depth  of  the  zone  of  P  enrichment  in  two  soil  pro¬ 
files,  Windsor  and  Charlton,  to  which  wastewater 
had  been  applied  for  5  years.  In  the  case  of  the 
Fort  Devens  soil,  field  measurements  of  P  enrichment 
were  in  agreement  with  the  prediction  that  the  P 
storage  capacity  of  this  profile  had  been  saturated 
and  that  wastewater  was  passing  the  163-cm  depth 
essentially  unchanged  with  respect  to  P  concentra¬ 
tions. 

Although  the  model  supplies  only  one  item  of 
information,  in  many  practical  situations  an  estimate 
of  the  P  storage  capacity  of  a  soil  profile  may  be  all 
that  is  required  when  evaluating  a  site  for  wastewater 
treatment.  The  model  can  be  used  to  evaluate  not 
only  proposed  treatment  areas  but  also  the  condition 
of  existing  treatment  areas,  providing  that  soil  ma¬ 
terials  representative  of  the  treatment  area  before 
wastewater  application  are  available.  The  latter  may 
generally  be  obtained  from  land  contiguous  with  the 
treatment  area. 

Although  the  generalized  P  sorption  isotherm  can 
be  used  to  predict  solution  P  concentrations,  the  P 
sorption  submodel  is  less  successful  in  the  prediction 
cf  soil  solution  P  concentration  at  a  particular  stage 
after  the  initiation  of  wastewater  treatment.  This 
probably  reflects  the  fact  that  the  P  sorption  sub¬ 


model  is  based  only  on  considerations  of  P  sorption 
and  does  not  include  the  impact  of  desorption  hys¬ 
teresis  (Barrow  and  Shaw  1975,  Ryden  and  Syers 
1977b)  in  controlling  soil  solution  P  concentration. 

The  model  also  assumes  a  sharp  boundary  between 
the  Z'  ?  of  P  enrichment  and  the  nonenriched  soil 
below.  Preferential  movement  of  solution  P  through 
macropores  and  channels  (Kanchanasut  et  al.  1978, 
Thomas  and  Phillips  1979)  within  the  profile  will 
result  in  a  more  rapid  movement  of  solution  P  to 
greater  depths  than  predicted  by  the  present  model. 

This  limitation  is  shared  with  other  models  in  that 
they  all  assume  uniform  movement  of  water  and  sol¬ 
ute  through  the  profile.  Nevertheless,  field  data  (Tables 
4  and  5)  obtained  for  low-rate  treatment  on  the  Wind¬ 
sor  and  Charlton  soils  suggest  that  the  assumption  of 
a  sharp  boundary  may  be  reasonably  valid. 

The  P  sorption  submodel  used  in  the  present  study 
assumes  that  sorption  is  the  only  process  involved 
in  P  retention  by  soils.  This  assumption  is  probably 
valid  for  a  wide  range  of  soil  types,  as  discussed  by 
Syers  and  Iskandar  (1980).  However,  on  calcareous 
soils  and  at  high  solution  P  concentrations,  it  is  pos¬ 
sible  that  precipitation  of  calcium  phosphates  may 
make  a  significant  contribution  to  P  retention  (Hol- 
ford  and  Mattingly  1975).  In  such  circumstances, 
data  may  not  conform  to  the  generalized  sorption 
isotherm.  Furthermore,  extraction  of  soils  with  NaOH 
does  not  remove  precipitated  calcium  phosphates, 
and  the  model  may  lead  to  an  underestimation  of  the 
P  storage  capacity. 

The  present  model  has  several  practical  advantages, 
a  major  one  being  that  the  P  sorption  submodel  re¬ 
quires  only  a  small  number  of  laboratory  measure¬ 
ments.  If  it  is  assumed  that  P  sorption  by  a  particular 
soil  conforms  to  the  generalized  P  sorption  isotherm, 
the  storage  capacity  of  a  single  depth  element  can  be 
determined  from  only  two  laboratory  measurements 
and  an  estimate  of  the  bulk  density  of  the  undisturbed 
soil.  As  an  increasing  number  of  soils  appear  to  con¬ 
form  to  the  generalized  P  sorption  isotherm,  it  may 
be  possible  in  future  applications  of  this  model  to 
eliminate  the  procedure  for  establishing  conformity 
to  the  generalized  isotherm  (Table  1).  The  assess¬ 
ment  of  P  sorption  at  equilibrium  using  extrapolation 
to  infinite  time  eliminates  the  impact  of  soil-to-solu- 
tion  ratio  on  the  estimate  of  P  sorption  by  soils  (Hope 
and  Syers  1976).  The  use  of  soil  suspensions  and  a 
wide  soii-to-solution  ratio  provides  a  simpler  method 
of  assessing  P  sorption  by  soils  than  the  use  of  sat¬ 
urated  soil  pastes,  soil  cores,  or  soil  columns. 

The  model  developed  in  the  present  study  offers 
potential  for  a  simple  estimation  of  the  P  storage 
capacity  of  the  soil  profile.  Initial  application  of 
the  model  has  produced  encouraging  results,  but 


further  evaluation  is  required  on  other  sites  before 
its  general  applicability  can  be  established. 

CONCLUSIONS 

The  P  model  evaluated  in  the  present  study  ade¬ 
quately  predicted  site  longevity  in  terms  of  P  re¬ 
moval  from  applied  wastewater.  The  model  has 
been  successful  in  predicting  the  lower  depth  of 
the  zone  of  P  enrichment  in  two  soil  profiles,  Wind¬ 
sor  and  Charlton,  to  which  wastewater  had  been 
applied  for  5  years.  In  the  case  of  the  Fort  Devens 
soil  (gravelly  sand),  field  measurements  of  P  enrich¬ 
ment  were  in  agreement  with  the  prediction  that 
the  P  storage  capacity  of  this  profile  had  been  sat¬ 
urated  and  that  wastewater  was  passing  the  163-cm 
depth  with  P  concentrations  essentially  unchanged. 
The  P  model  is  particularly  useful  in  selecting  and 
evaluating  sites  for  land  treatment,  as  well  as  for 
site  management. 

RESEARCH  NEEDS 

As  indicated  above,  the  model  evaluated  in  the 
present  study  has  previously  received  only  limited 
evaluation.  The  P  sorption  submodel  requires  more 
extensive  evaluation  using  a  wide  range  of  soil  types 
from  different  localities.  The  applicability  of  the 
submodel  will  be  most  easily  established  by  testing 
equilibrium  P  sorption  data  for  conformity  with  the 
generalized  P  sorption  isotherm.  In  particular,  the 
P  sorption  submodel  requires  further  evaluation  using 
calcareous  soils  to  evaluate  the  possible  importance 
of  precipitation  in  P  removal  by  these  soils.  The  use 
of  the  P  sorption  submodel  to  predict  storage  capac¬ 
ity  and  the  P  balance  model  to  predict  site  longevity 
requires  use  and  evaluation  at  a  wide  range  of  waste- 
water  treatment  sites.  This  evaluation  will  be  most 
valuable  if  conducted  at  treatment  sites  that  have 
been  in  operation  for  at  least  5  to  10  years  and  where 
soil  samples,  representative  of  the  profile  wastewater 
addition,  are  available.  Information  obtained  in 
these  studies  will  define  the  scope  of  applicability  of 
this  simple  model  to  determine  storage  capacity  and 
site  longevity  with  respect  to  P  removal  from  waste- 
water. 
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